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Abstract 
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fractional  combinatorial  optimization.  We  prove  that  Newton’s  method  runs 
in  0(m)  iterations  for  the  maximum  mean- weight  cut  problem.  One  iteration 
is  dominated  by  the  TnaYiTnuTn  flow  computation,  so  the  overall  running  time 
is. O(m^n).  The  previous  fastest  algorithm  is  bsised  on  Meggido’s  parametric 
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1  Introduction 


A  linear  fractional  combinatorial  optimization  (LFCO)  problem  is  defined  as  foHows.  A 
specification  of  a  system  X  C  {0,  Ip  and  two  vectors  (oi,  02, Op),  (6i,  62,  •••»  ip)  G 


za'e  given. 


MaTimize 


subject  to  (*i, *2}  •  •  •  5 ®p)  €  X. 


a-iXi  -|-  02^2  +  •  •  •  +  Qp^p 

tiajj  +  ^>2*2  +  ■  — I"  ip®p 
X  belongs  to  some  specified  class  of  systems.  Vector  (®i,a:2,  ...,Xp)eX,  represents  some 
structure,  a  subset  of  some  underlying  set  of  p  objects,  a,  and  bi  are  the  cost  and  the 


weight  of  object  i.  YLi  (SLi  ^iXi)  axe  the  cost,  weight, 

and  mean-weight  cost  of  the  structure  represented  by  (si,  *2)  •  •  •  >  *p)-  Ill  an  LFCO  problem 
we  are  asked  to  compute  the  maximum  mean-weight  cost  of  a  structure  in  X .  We  also  want 
to  find  a  structure  which  has  the  maximum  mean-weight  cost. 


For  example  the  TnaTiTnuTn  mean-weight  spanning  tree  problem  is  an  LFCO  problem. 
In  this  problem  the  class  of  systems  is  the  class  of  sets  of  all  spanning  trees  in  graphs: 

{{T I T  —  spanning  tree  iaG}  \  G  —  graph}. 

The  description  of  graph  G  =  (V,F)  is  the  specification  of  Xq  Q  {0,  l}*",  the  set  of  (the 
characteristic  vectors  of)  all  spanning  trees  in  G  (m  =  |F?1).  The  underlying  set  of  objects 
for  Xg  is  the  set  of  edges  in  G.  Vector  (xi,  xj, . . . ,  x„)  €  Xa  represents  spanning  tree  T  in 
G  such  that  edge  i  belongs  to  T  iff  Xf  =  1.  Numbers  ai,  02, . . . ,  Om  axe  the  costs  of  edges  and 
61, 62,  •  •  • ,  axe  the  weights.  The  task  is  to  find  a  spanning  tree  in  G  with  the  maximum 
mean- weight  cost. 

We  will  use  bold  letters  to  denote  vectors  and  corresponding  italic  letters  to  denote 
their  coordinates,  for  example  x  =  (xj ,  X2,  •  •  • ,  ®p)-  For  two  vectors  a  and  x  of  equal  length, 

ax  denotes  the  scalar  product  of  a  and  x,  t.c.,  ax  =  oiXi  -I-  02X2  -1 - 1-  ap®p*  With  this 

notation  an  LFCO  problem  is  to 

ax 

maximize  subject  to  x  E  Af. 

We  assume  that  bx  >  0  for  all  x  €  -V,  and  ax  >  0  for  some  x  €  -V.  Problem  (PI)  can  be 
equivalently  formulated  in  the  following  way: 

(P2)  minimize  6^  subject  to  (a^c)  —  5(bx)  <  0,  for  all  x  £  Af. 

Let  6*  denote  the  miTiiTrmTn  we  are  looking  for.  Let 


Function.  h{S)  is  convex,  piecewise  linear,  decreasing,  and  S’  is  its  only  root.  We  assume 
that  we  can  maximize  any  linear  function  over  X,  so  for  2iny  6  we  can  compute  .h{S)  and 
X  €  A’  such  that  (a  —  ^b)x  =  h{S).  Let  >1  be  an  algorithm  which  does  such  computation 
zind  let  T  be  its  running  time.  Having  such  an  algorithm  enables  us  to  compute  S"  with 
Newton’s  method.  We  show  that  Newton’s  method  finds  6“  in  O(p^)  iterations.  This 
bound  on  the  number  of  iterations  is  independent  of  the  complexity  of  the  cleiss  of  systems 
which  X  belongs  to,  and  independent  of  the  complexity  of  computing  h{6).  One  iteration 
is  dominated  by  algorithm  A,  so  the  overall  time  is  0{p^T). 

The  other  way  of  computing  S'  is  Megiddo’s  parametric  search  method  [8].  This  ap¬ 
proach  gives  an  algorithm  with  running  time  0{TQ{T'logQ  -1-  Q)),  where  T'  and  Tq  refer 
to  the  sequential  running  time  and  the  parallel  running  time  on  Q  processors  of  computing 
the  sign  of  h{6)  for  a  given  6.  Thus  an  efficient  algorithm  for  an  LFCO  problem  can  be 
obtained  if  an  efficient  parallel  algorithm  for  the  corresponding  nonfractional  problem  is 
available. 

If  costs  and  weights  are  integers  firom  [—17,  U],  then  using  the  straightforward  binary 
search  we  can  find  in  0{log{pU))  iterations  an  interval  I  of  length  less  than  lf{pU)^  which 
contains  S'.  One  iteration  is  dominated  by  computing  the  sign  of  k{S)  for  2ui  appropriate 
S.  An  interval  of  length  less  than  l/{pU)^  can  contain  at  most  one  mean-weight  cost.  The 
overall  time  of  this  approach  is  0(r'log(pl7)  -f  T). 

In  the  second  part  of  the  paper  we  consider  the  minimum  maximum  arc  cost  trans¬ 
shipment  (MAC)  piroblem,  which  is  the  transshipment  problem  when  the  maximum  of  the 
arc  costs  is  being  minimized.  This  problem  can  be  reduced  to  the  maximum  mean-weight 
cut  (MWC)  problem,  which  is  an  LFCO  problem.  The  underlying  nonfractional  problem 
is  the  standard  maximum  flow  problem.  This  problem  does  not  have  an  efficient  parallel 
algorithm.  Nevertheless,  up  to  now,  Megiddo’s  method  gave  the  best  known  upper  botmd, 
namely  6(mn®)  [8],  for  the  MWC  and  MAC  problems.  We  show  that  Newton’s  method 
applied  to  the  MWC  problem  runs  in  0(m)  iterations.  One  iteration  is  dominated  by  max¬ 
imum  flow  computation,  so  the  overall  running  time  is  O(Tn^n).  The  same  boimd  holds  for 
the  MAC  problem. 

2  Newton’s  Method  for  Linear  Fractional  Combinatorial 
Optimization 

The  following  is  Newton’s  method  for  computing  the  root  S'  of  h{S)  (defined  in  (1)). 
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Figure  1:  Newton’s  method  for  solving  —  0 


Let  ^  be  the  current  estimation  to  6".  InitiaUy  6  =  0.  (Observe  that  our  as¬ 

sumptions  guarantee  S'  >  0.)  During  one  iteration  we  compute  fc(^)  and  i  €  such  that 
h(S)  =  (ax)  -  ^(bx),  t.c.,  we  maximke  linear  function  (a  -  ^)x  over  X.  This  is  done  by 
algorithm  A  If  h{6)  =  0,  then  S'  =  6  and  the  algorithm  taminates.  Otherwise  we  compute 
the  next  estimation  6  *-  ax/bx,  the  mean- weight  cost  of  x,  and  go  to  the  next  iteration. 
The  process  is  illustrated  in  Figure  1. 

Let  Si  be  the  value  of  6  at  the  beginning  of  ith  iteration,  and  Xi,  Hi,  and  Bi,  be  x, 
(a  —  ^fb)x,  and  bx  from  this  iteration. 


Hi 

Hi 


(a  -  tfib)xi 


bxi, 

axi 


bxi’ 


m£a:{(a  —  5ib)x  |  x  €  X}, 


and  it  can  be  easily  derived  that 


(2) 
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The  following  lemma  indicates  fast  convergence  of  the  above  algorithm. 


Lemma  1 


Ei  ^  Bi  -  ' 


Proof.  Vector  Xi  maximizes  (a  -  ^,b)x,  so 


(a  -  5fb)xi  >  (a  -  ^ib)Xt+i. 


Therefore 


Ei  =  (a  -  ^,b)xi  >  (a  -  5,b)xi+i 

=  (a  —  ^t+ib)xi4.i  +  (5»+i  —  ^j)bxf+i 

=  Ei+i  + 

This  implies  Inequality  (3).  ■ 


The  above  lemma  gives  immediately  the  following  bound  on  the  number  of  iterations. 

Theorem  1  If  in  an  LFCO  problem  the  coordinates  of  vectors  a  and  b  are  integers  from 
[-17,  U],  then  Newton’s  method  finds  the  optimum  in  0(log(pl7))  iterations. 


Proof.  Inequality  (3)  implies 


For  any  i  for  which  Ei  >  0, 


Ei^\Bi+i  ^  1 

EiBi  -  4‘ 


{pUf  >  EiBi  > 


pU' 


The  second  inequality  holds  because 


_  ax.bx._i  -  axf-ibxt  J_ 

'  bxi_i  -pU' 


Inequalities  (4)  and  (5)  imply  the  bound  of  0(log(pCf))  on  the  number  of  iteration. 


(4) 

(5) 


The  special  case  with  uniform  weights,  t.c.,  b  =  (1, 1, . . . ,  1),  is  of  independent  interest. 
The  task  here  is  to  find  the  maximum  mean  cost  of  a  structure.  Sequence  (Bi),  excluding  the 
last  iteration,  is  decresising.  This  implies  the  following  strongly  polynomial  bound,  which 
was  observed  by  Karzanov  [11]  and,  in  the  context  of  the  maximum  mean  cut  problem,  by 
McCormick  and  Ervolina  [7]. 

Theorem  2  If  the  weights  are  uniform,  then  Newton’s  method  runs  in  at  most  p+1  iterations. 
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Now  we  prove  a  strongly  polynomial  bound  on  the  number  of  iterations  for  the  general 
case,  when  both  costs  and  weights  can  be  arbitrary  real  numbers.  The  intuition  behind 
the  analysis  is  as  follows.  Lemma  1  suggests  that  there  are  some  sequences  related  to  the 
rate  of  convergence  of  the  method,  which  tend  to  zero  at  least  geometrically  fast.  The 
elements  of  these  sequences  are  obteiined  from  only  2p  numbers,  the  costs  ai,a2, . .  .,0^,, 
and  the  weights  bi,b2..-,bp,  using  only  0(p)  additions/subtractions  and  at  most  one  or 
two  multiplications/divisions.  We  show  that  because  of  this  limiting  use  of  arithmetical 
operations,  the  lengths  of  such  sequences  cannot  be  too  big. 

To  expand  this  intuition  a  little  bit,  let  us  assume  that  Bi+i/Bi  <  1  -  a,  for  some 
positive  constant  a,  that  is  sequence  (5,)  tends  to  zero  at  least  geometrically  fast.  Each 
element  of  sequence  (5i)  is  a  sum  of  different  elements  from  {61,  b2,...,  bp}.  Let  us  further 
assume  that  bi  >  b2  >••>  bp  >  0.  Obviously  Bi  <  pbi.  Since  (£»)  decreases  at  least 
geometrically,  for  some  I  =  O(logp),  Bi  <  b^.  It  means  that  61  is  not  a  term  in  Bi,  nor  is 
it  in  Bi,  for  any  i  >  1.  Bi  <  {p  -  1)^2?  and  the  next  O(logp)  elements  of  (B,)  exclude  62. 
then  ba,  and  so  on.  Therefore  the  length  of  sequence  (Bi)  is  O(plogp). 

There  are  two  reasons  why  the  general  ceise  is  more  complicated.  First,  we  have  to  deal 
with  positive  and  negative  nTimbers.  Even  if  both  costs  and  weights  are  positive,  negative 
numbers  appear  because  subtractions  are  used  in  forming  the  elements  of  (Bi).  Second,  if 
(Bi)  does  not  decrease  fast  enough,  then  we  have  to  use  sequence  (Bi),  whose  elements  are 
not  just  subsums  of  some  given  set  of  numbers.  The  following  lemma  is  our  tool  to  deal 
with  both  positive  and  negative  numbers.  Here  the  coordinates  of  vector  c  are  numbers 
which  are  used  to  form  sums,  and  (ytc)  is  a  sequence  of  such  sums. 

Lemma  2  Let  c  6  and  €  {0, for  fc  =  1, 2, . . . ,  g,  be  such  that 

|yfc+ic|  <  ^^ly^c],  for  i;  =  1, 2, . . ., g  —  1.  (6) 

Then  q  <p. 

Proof.  Assume  that  q  >  p.  Let  Y  denote  the  matrix  whose  rows  are  vectors  yi,  y2,  • . ., 
Yg.  Since  q  >  p,  there  is  a  nonzero  vector  z  €  R®  such  that  zY  =  0.  Moreover,  z  can  be 
chosen  in  sudi  a  way  that  all  its  coordinates  are  rational  numbers  with  denominators  and 
enumerators  being  the  determinants  of  submatrices  of  Y.  Let  z  =  (21,  Z2»  •  •  •  >  ^g)  such  a 
vector.  Thus,  for  t  =  1, 2, . . . ,  g,  either  Zf  =  0  or 
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We  have 


0  =  zYc  =  2l(yic)  +  22(y2C)  + - h  2g(y9C). 

But  this  is  not  possible,  because  the  nonzero  terms  on  the  right  side  cannot  cancel.  Since 
y,  7^  0,  so  for  some  1  <  i  <  ff  -  1,  z,-  7^  0,  and  2iy,-c  #0.  E  1  <  i  <  j  <  g  and  Ziy.-c  7^  0, 
then  (6)  and  (7)  imply 

|2i(yic)l  =  ^k(y.-c)|. 


Lemma  3  There  are  at  most  O(p^logp)  iterations  k  such  that  Bk+i  <  \Bk- 

Proof.  Consider  the  sequence  of  all  iterations  k  such  that  <  fPfc-  Take  from  this 

sequence  every  /th  iteration,  where  I  =  fbplogp].  Let  they  be  iterations  •  •  •  ,iq-  For 

1  <  *  <  ?  -  1, 

<  (3)  -Fit  < 

Lemma  2  implies  g  <  p,  which  implies  0{p^  logp)  bound  on  the  number  of  iterations  k  such 
that  Bk+i  <  iBk.  m 

Lemma  4  There  are  at  most  C?(p^logp)  consecutive  iterations  k  such  that  Bk+i  >  §£*• 

Proof.  Consider  a  sequence  of  consecutive  iterations  t,  t  +  1,  . . j,  such  that  for  each 
i<  k  <  j  —  1: 

Bk+i  >  ^Bk.  (8) 

This  and  Inequality  (3)  imply  that  for  each  i  <  k  <  j  —  1: 

Bk+i  < 


Therefore 


-  fM  =  ^  ~  -  &)■ 


(9) 


Let  I  =  [3plogp]  +  1.  Consider  the  sequence  of  every  /th  iteration  out  of  iterations  t, 
t  +  1,  . . j.  Let  they  be  iterations  *1,  *29  •  •  •  >  *?•  We  show  that  Lemma  2  can  be  applied 
to  the  sequence  (xi4_i(-a  +  where  1  <  Ife  <  g.  This  is  a  sequence  of  subsums  of 

a,  +  Si^br).  First  we  estimate  the  rate  of  convergence  of  sequence  {Si^)  to  6i^.  For 

l<jb<g-2, 


<  (^i,  -  -  ^1,-2)  H - h  (^it+i+l  - 
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<  (^i*+i+i  -  ^  i  +  ■  ■  ■  +  ) 

<  2(5ij^j+i  —  ^t*+i) 

/T  \  (ijk+x-it) 

£  2  J  (^it+l  ~  ^»*) 

The  inequalities  in  the  second  and  fourth  lines  follow  from  Inequality  (9). 
Recall  that  6i+i  =  {axi)/{hxi),  so  for  t  =  1, 2, . . . , 

Xi,_a(-a+  5i,b)  =  (5i,  - 


Put 

c  =  -a  +  5i,b, 

yk  =  for  fc  = 

We  have  for  t  =  1, 2, . . . ,  g  —  1, 

yife+ic  =  -Bit+j -1  (^i,  -  ^it+, )  <  -Bit-i ^ (^i,  -  )  =  ^y*®- 

Thus  vectors  c  and  yjt’s  satisfy  the  conditions  of  Lemma  2,  so  g  <  p.  It  means  that  there 
are  at  most  O(p^logp)  consecutive  iterations  for  which  Inequality  (8)  holds.  ■ 

Lemmas  3  and  4  immediately  imply  the  following  strongly  polynomial  bound  on  the 
number  of  iterations. 

Theorem  3  Newton's  method  applied  to  an  LFCO  problem  finds  the  optimum  in  O(p^log^p) 
iterations. 

If  the  weights  are  nonnegative,  then  there  are  at  most  O(plogp)  iterations  with  < 
Therefore  in  such  a  case  Newton’s  method  works  in  0{p^)  iterations. 

3  Minimum  maximum  arc  cost  transshipment  problem 


A  (transshipment)  network  G  =  (V,  £,«,</)  is  a  digraph  (V^)  with  a  capacity  function 
u  ;  E  - — *■  R  U  {+oo},  and  a  demand  function  d  :  V  — »  R  such  that  53«€V 
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assume,  for  convenience  in  presentation  but  without  losing  generality,  that  E  is  symmetric 
and  u  is  nonnegative.  Negative  d(v)  means  that  w  is  a  source  -  a  vertex  with  supply,  positive 
d{v)  means  that  v  is  a  sink  -  a  vertex  with  demand.  Let  n  and  m  denote  the  cardinality  of 

V  and  the  half  of  the  cardinality  of  E,  respectively.  We  assume  m>n. 

We  adopt  the  convention  that  any  function  ip  :  F  — *  R  is  extended  to  all  finite  subsets 
of  F:  (p{A)  =  fp{x).  Thus  if  A  C  V,  then  d{A)  =  d[v),  which  is  the  net  demand 

in  A.  The  total  demand  in  G  is  d{G)  =  J^v:d{v)>o 

A  pseudoflow  in  G  is  an  antisymmetric  function  f  :  E  — >  R,  i.e.,  for  every  (r,  w)  £  E, 
f{v,w)  =  —f{w,v).  A  pseudoflow  /  creates  excesses  at  vertices;  the  excess  at  vertex  v  is 
el{v)  =  Yiwiwv)€E  ^  “  ^(®))  pseudoflow  f  satisfies  the 

demand  at  vertex  v.  A  flow  is  a  pseudoflow  f  that  meets  capacity  constraints  and  does  not 
increase  demands:  f{w,v)  <  u(t[;,t;),  for  every  arc  (w,v)  €  E,  and  e^(v)  is  between  0  and 
d(v),  for  every  v  £  V.  A  transshipment  is  a  flow  which  satisfies  demands  at  all  vertices. 
G  is  feasible  if  there  is  a  transshipment  in  G.  The  transshipment  feasibility  problem  is  the 
problem  of  verifying  if  a  given  network  G  is  feasible.  For  a  pseudofiow  f  in  G,  the  residual 
network  is  {V,E,u^  =  u  —  f,d^  =  d—  e^).  Ga  denotes  G  with  the  capacity  function 
changed  to  u.  Thus  Gg  =  {V,E,u^,d^). 

If  5  and  T  partition  V,  then  cut  {S,T)  in  G  is  the  set  of  arcs  {v,w)  of  E  such  that 

V  £  S  and  w  £T.  The  capacity  and  the  surplus  of  a  cut  {S,T)  are  u{S,T)  =  2e€(5,r)“(®) 
and  surpluses, T)  =  d{T)  —  u{S,T),  respectively.  Our  definition  of  surplus{S,T)  is  the 
same  as  in  [7]  the  definition  of  V{T,  S),  the  v^due  of  a  cut  (T,  S).  It  is  easy  to  verify  that 
the  surplus  of  a  cut  is  the  same  in  G  as  in  G^  for  any  flow  /.  A  positive  surplus  means 
that  the  cut  blocks  the  flow.  Any  flow  must  leave  at  5  at  least  the  amount  surplus{S,T) 
of  the  commodity,  which  is  demanded  at  T.  Let  mean{S,T)  =  surp/tts(5,r)/|(5,r)| 
be  the  mean  surplus  of  (5,  T)  in  G,  2md,  if  a  weight  function  b  :  E  — ►  R  is  specified, 
meanjw{S,T)  =  surplus{S,T)/b{S,T)  be  the  mean-weight  surplus  of  (5,T)  in  G. 

The  minimum  cost  transshipment  problem,  which  is  often  refer  to  as  the  transshipment 
problem,  is  to  compute  for  a  given  network  G  and  a  given  cost  function  c  :  E  — *  R,  the 
mmiTnum  cost  of  a  transshipment  in  G,  when  the  cost  of  transshipment  /  is 
This  problem  has  been  well  studied,  often  as  the  minimum  cost  circulation  problem  (these 
problems  are  reducible  to  each  other  in  linear  time).  Fastest  currently  known  algorithms  for 
this  problem  can  be  found  in  [1,  2,  4,  9, 12].  In  this  paper  we  are  interested  in  the  problem 
of  minimizing  the  TnaTiTnum  arc  cost  instead  of  the  sum  of  arc  costs.  Here  is  the  list  of  the 
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transshipment  problems  we  consider.  They  are  closely  related  to  each  other.  We  assume 
that  cost  functions  and  weight  functions  are  noimegative. 


Minimum  maximum  arc  cost  transshipment  problem  (MAC) 

A  transshipment  network  G  and  a  cost  function  c  :  E  — >  R  are  given.  Compute  the 
TnTnirrmTn  of  the  m a-TiTmim  3xc  cost  over  all  transshipments  in  G,  that  is,  compute 

6"  =  inin{max{/(e)c(e)}  |  /  is  a  transshipment  inG}. 

Minimum  capacity  violation  cost  transshipment  problem  (MCVC) 

A  transshipment  network  G  is  given,  which  might  not  be  feasible.  The  objective  is  to 
obtain  a  feasible  network  by  increasing  the  capacities.  A  weight  function  b  :  E  — >  R 
is  given,  which  specifies  by  how  much  the  capacities  of  the  arcs  can  be  increased  in 
unit  cost.  It  means  that  the  cost  of  increase  of  the  capacity  of  axe  e  from  «(e)  to 
f  is  i/b{e).  Compute  the  Tni-ni-rmiTn  of  the  maximum  cost  of  increcising  the 
capacity  of  an  arc,  th?it  is,  compute 

S"  =  min{^  >  0  |  feasible}. 


Minimum  maximum  capacity  violation  treinsshipment  problem  (MCV) 

This  is  the  MCVC  problem  with  the  uniform  weight  function:  6  =  1. 

Mcocimum  mean- weight  cut  problem  (MWC) 

A  transshipment  network  G  and  a  weight  function  b  :  E  — R  are  given.  Compute 
the  maximum  mean-weight  surplus  over  all  cuts  with  noimegative  surplus,  that  is, 
compute 

S"  =  max{  I  {S,  T)  cut  in  G  with  nonnegative  surplus}. 

6(5,  T) 


Maximum  mean  cut  problem  (MC) 

This  is  the  MWC  problem  with  the  uniform  weight  function.  It  means  that  we  are 
asked  to  find 


6*  =  max{ 


surplus{S,  T) 

1(5,  r)| 


1  (5,  r)  cut  in  G  with  nonnegative  surplus}. 


A  cut  with  the  TnaTimTiTn  nonnegative  mean-weight  surplus  and  a  cut  with  the  maxi¬ 
mum  nonnegative  mean  surplus  we  will  call  a  Ttiaximum  mean-weight  (surplus)  cut  and  a 
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maximum  mean  (surplus)  cut,  respectively.  We  formulated  the  above  problems  in  a  way 
that  the  TniTiirmiTn  or  the  TnaTiTmiTu  veilue  is  sought,  not  the  corresponding  optimal  flow  or 
cut.  We  did  it  only  for  convenience  in  farther  presentation.  The  algorithm  we  discuss  in 
the  next  section  and  actually  all  other  algorithms  for  the  above  problems  compute  both  the 
optimal  value  and  the  corresponding  optimal  structure  (transshipment  or  cut). 

As  mentioned  in  [8],  the  problem  of  minimizing  the  maximum  arc  cost  appears,  for 
example,  when  it  is  desirable  not  to  use  an  arc  at  its  full  capacity.  Each  arc  e  has  some 
safety  level  s(c)  and  the  ratio  f{e)ls{e)  is  sought  to  be  minimized,  with  the  general  objective 
of  getting  the  maximal  ratio  as  small  as  possible.  The  MAC  problem  captures  also  the 
following  dynamic  situation.  The  network  is  in  continuous  activity,  that  is,  the  sources 
continuously  keep  producing  commodity  and  sending  it  to  the  sinks.  The  objective  is  to 
maximize  the  throughput.  More  formally,  suppose  l/c(e)  is  the  amount  of  the  commodity 
which  can  be  shipped  through  arc  e  in  unit  time.  It  means  that  the  shipment  of  /(e)  units 
through  e  takes  /(e)c(e)  time.  Further  suppose  that  every  source  s  generates  — d(s)  units 
of  the  commodity  every  t  units  of  time  and  sends  it  out,  and  every  sink  t  has  to  get  d(t) 
units  of  the  commodity  every  r  units  of  time.  Then  the  optimal  6’  in  the  formulation  of 
the  MAC  problem  is  equal  to  the  minimal  possible  t  such  that  no  congestion  arises  at  any 
vertex. 

Before  discussing  the  previous  work  relevant  to  the  problems  in  which  we  axe  interested, 
we  first  establish  relation  between  them. 

The  value  of  a  flow  f  is  the  amount  by  which  /  decreases  the  total  demand:  value{f)  = 
I3v:d(tr)>o  (®)-  maximum  flow  in  G  is  a  flow  in  G  of  the  maximum  value.  The  maximum 
flow  problem  is  the  problem  of  finding  a  mstximum  flow  in  a  given  transshipment  network. 
A  maximum  (surplus)  cut  is  a  cut  with  the  maximum  noimegative  surplus.  The  maximum 
(surplus)  cut  problem  is  the  problem  of  finding  a  maximum  cut  in  a  given  transshipment 
network. 

The  transshipment  feasibility  problem,  the  maximum  flow  problem  and  the  maximum 
cut  problem  are  the  non&actional/nonparametric  problems  that  correspond  to  our  frac¬ 
tional/parametric  problems.  The  transshipment  feasibility  problem  and  the  maximum  flow 
problem  can  be  easily  reduced  to  the  standard  maximtim  flow  problem.  A  maximum  flow 
in  a  transshipment  network  can  be  used  to  find  a  maximum  surplus  cut,  in  a  similar  way  as 
a  TTiaTiTniTm  flow  in  a  standard  TnaviTnTim  flow  network  can  be  used  to  identify  a  minimum 
capacity  cut.  Let  TMF(n,m)  denote  G(M(n,  m)),  where  M{n,m)  is  the  time  complex- 
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ity  of  the  standard  maximum  flow  problem.  The  transshipment  feasibility,  the  maximum 
flow  in  a  transshipment  network,  and  the  maximum  surplus  cut  problems  can  be  solved  in 
TMp{n,  m)  time.  Currently  it  is  known  that  TMF{n,m)  =  O(nmlog(n^/m))  [3],  and  slightly 
better  bounds  exist  for  dense  graphs. 

Let  for  a  transshipment  problem  P,  rp(n,  m)  denote  its  time  complexity.  The  following 
theorem  shows  close  relation  between  the  MAC  problem  and  the  MCVC  problem. 

Theorem  4 

1-  ?MCVc(^)  ^  2m). 

2.  TuAc{.n,m)  <  O(mlogm)  +  TMF(n,»n)0(logm)  +  Tucvc{n,rn). 

Proof. 

1.  For  a  given  instance  I  of  the  MCVC  problem  construct  an  instance  I'  of  the  MAC 
problem  in  the  following  way.  Replace  each  arc  e  in  I  with  two  parallel  arcs  ej  and 
C2.  Put  «'(ei)  =  u(c),  c'(ei)  =  0,  tt'(c2)  =  +oo,  c'(e2)  =  l/6(e).  The  demand  function 
in  r  is  the  same  as  in  J.  A  solution  (transshipment)  to  I'  gives  in  the  natural  way 
the  solution  to  I. 

2.  Let  be  the  solution  to  the  given  instance  I  of  the  MAC  problem.  It  is  easy  to  see 

that  r  <  ^  if  and  only  if  «  feasible.  Therefore 

S"  =  min{5  >  0  |  C?nun{u,«/c}  **  feasible}. 

To  get  an  instance  of  the  MCVC  problem,  it  is  enough  to  know  for  each  arc  c,  if 
u(e)  is  essential,  i.e.,  if  «(c)  <  ^’/c(c).  We  sort  the  list  {u(c)c(e)  \ee  E}  and  using 
binary  search  we  find  the  position  of  S*  in  it.  This  takes  Turin,  m)0(\ogm)  time. 
Now  we  construct  an  instance  T  of  the  hfCVC  problem.  The  underl3dng  graph  and 
the  demand  function  remain  without  changes.  The  capacities  and  weights  are  defined 
as  follows. 

iv!{e),h\e))  = 

6*,  the  solution  to  I,  is  also  the  solution  to  V. 


/  (u(c),0)  if  u(c)c(c)  <  r, 
\  (0,  l/c(c))  if  tt(e)c(e)  >  T. 
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Let  meanjws{S,T)  denote  the  mean-weight  surplus  of  (5,  T)  in  Then 

meanjws{S,T)  =  meanjw{S,T)  —  6. 

It  means  that  the  maximum  mean_u?(5,  T)  is  equal  to  the  minimum  S  for  which  Gt,^sb  does 
not  contEun  cuts  with  positive  surplus.  Therefore  the  following  theorem  imphes  that  the 
MCVC  problem  and  the  MWC  problem  are  equivalent  (and  so  are  MCV  and  MC),  i.e., 

min{5  >  0  I  Gu+ifc  is  feasible}  = 

m.ax{mean.w{S,  T)  |  (5,  T)  cut  in  G  with  nonnegative  surplus}. 

Theorem  5  [5]  G  is  feasible  iff  there  is  no  cut  in  G  with  positive  surplus. 

The  MWC  problem  is  an  LFCO  problem.  Xg  C  {0, 1}"+™  corresponds  to  the  set  of  all 
cuts  in  G.  Vector  x  €  Xg  represents  cut  {S,T)  such  that  for  1  <  i  <  n,  *,•  =  1  iff  Wj  G  T, 
and  for  I  <j  <m,  Xn+j  =  1  iff  €  (5,  T),  where  Vi  is  the  ith  vertex  and  Cj  is  the  jth  arc. 
Vector  a  from  the  formulation  of  an  LFCO  problem  is  here  (d,  -u),  where  d  and  u  are  the 
vector  representations  of  functions  d  and  u.  If  x  represents  (5,  T),  then  ax  and  bx  are  now 
surpluses,  T)  and  6(5,  T). 

Megiddo  [8]  showed  how  to  solve  the  MAC  problem  with  his  parcunetric  search  method. 
This  method  gives  an  0(n®m)  running  time  and  can  give  a  better  time  only  if  a  break- 
throu^  in  parallel  maximum  flow  computation  is  achieved.  Binary  search  solves  MAC  with 
integral  data  in  0(TMF(«j»”)log(«t^))  time.  More  is  known  about  the  ca«e  with  uniform 
weights.  McCormidc  and  Ervolina  [7]  proposed  an  algorithm  for  computing  a  maximum 
TTieaTi  cut.  This  algorithm  is  essentially  an  application  of  Newton’s  method.  They  showed 
an  0{rn)  bound  on  the  number  of  iterations  (see  Theorem  2).  One  iteration  is  dominated 
by  Tna.TnmnTn  flow  computation.  We  were  able  to  prove  that  the  number  of  iterations  in 
this  algorithm  is  in  fact  0(n)  [10],  so  the  time  complexity  of  MC,  as  well  as  MCV  and 
MAC  with  uniform  costs,  is  0(n*m).  For  these  problems  the  bound  0(TMF(»i>^)los(^^)) 
can  be  improved  to  0{jim\og{nU))  by  using  approximate  maximum  flow  computation  in 
connection  with  either  Newton’s  method  [10]  or  binary  search  [6]. 

In  the  next  section  we  prove  an  0(w,)  bound  on  the  number  of  iterations  in  Newton’s 
method  applied  to  the  MCVC  problem.  This  implies  an  0(rMF(^j^)’’*)  =  0(m^n)  bound 
on  the  time  complexity  of  MCVC,  MAC  and  MWC. 


12 


4  Newton’s  method  for  the  MCVC  problem 


A  transshipment  network  G  and  a  (nonnegative)  weight  function  b  are  given.  We  assume 
that  G  is  not  feasible.  Ge  stands  for  and  generally  subscript  8  indicates  that  the 

underlying  network  is  G^+Bh-  We  want  to  find 

8’  =  min{^  |  G; is  feasible}, 

which,  as  mentioned  in  the  previous  section,  is  equal  to 

max{mean_uj(5,  T)  |  (5,  T)  cut  in  G  }. 

The  function  h{8)  used  in  the  discussion  of  a  general  LFCO  problem  in  Section  1  is  now 
h{6)  —  mzs.{surpluss{S,  T)  |  (5,  T)  cut  in  G  }. 

h(5)  is  also  equal  to  the  total  demand  left  in  Gg,  where  /  is  a  maximum  flow  in  G«. 

We  describe  the  way  Newton’s  method  finds  6“,  the  root  of  h{8).  Let  6  <  8*  he  the 
current  (under)estimation  of  8*.  Initially  6  =  0.  During  one  iteration  we  find  a  maximum 
cut  (5,  T)  in  G^.  This  is  done  by  computing  a  maximum  flow  /  in  G|.  If  the  surplus  of 
(5,  T)  is  zero,  i.e.,  if  /  is  a  transshipment,  then  8"  =  8  and  algorithm  terminates.  Otherwise 
we  compute  the  next  estimation  8  =  Tneanjw^^S,  T)  and  go  to  the  next  iteration.  Moving  to 
the  next  estimation  should  be  seen  as  increasing  the  capacity  so  that  the  surplus  of  {S,T) 
decreases  to  zero. 

For  any  flow  f  in.  G,  the  surpluses  of  cuts  are  the  same  in  G  as  in  G^ .  It  means  that 
in  the  current  iteration  we  can  perform  computation  on  any  residual  network  of  Gj  instead 
of  on  Gi  itself.  In  the  analysis  we  assume  that  the  network  in  the  current  iteration  is  G^, 
where  /  is  the  sum  of  the  TuaTiTmiTn  flows  computed  in  all  previous  itwations.  It  means 
that  /  is  a  maximum  flow  in  G;»,  where  S'  is  6  from  the  previous  iteration. 

Let  Si  be  S  at  the  beginning  of  iteration  i,  and  {Si,  Ti)  be  the  cut  found  in  this  iteration. 
Bi,  and  Si  from  Section  2  have  here  the  following  meaning. 

Bi  =  b{Si,Ti), 

Hi  =  h{Si)  =  surplus{Si,Ti)  -  Sib{Si,Ti). 

Let  fi  be  the  sum  of  the  flows  computed  through  iteration  i.  Flow  fi  is  a  maximum  flow  in 
Gsi-  Hi  is  the  total  demand  in  G^S 
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Theorem  1  implies  that  if  the  capacities  and  the  weights  are  integers  from  [—U,U], 
then  the  algorithm  runs  in  0(log(nCf))  iterations.  The  weights  are  nonnegative,  so  the 
remark  after  Theorem  3  implies  the  strongly  polynomicil  bound  of  (5(Tn^)  on  the  number  of 
iterations.  In  this  section  we  show  an  0{Tn)  bound  on  the  number  of  iterations. 

An  arc  e  is  unessential  in  G,  if  its  capacity  is  greater  than  the  total  demand  or  its 
weight  is  greater  than  the  weight  of  a  maximum  surplus  cut.  When  we  use  expression  “c  is 
unessential”  in  the  context  of  the  algorithm,  we  meein  that  e  is  imessential  in  the  current 
network.  U  an  arc  is  unessential  in  Gg,  for  some  6  >  and  flow  f  in  Gs,  then  it  is  unessential 
in  Gg^g! ,  for  any  6  >  0  and  any  flow  f  in  Gg^g,:  Increasing  6  increases  the  capacity  of 
an  arc  and  decreases  the  TnaYiTrmm  surplus  of  a  cut;  augmenting  with  a  flow  may  decrease 
the  capacity  of  an  arc  but  by  not  more  than  it  decreases  the  toted  demand.  Therefore  if 
at  some  point  in  the  algorithm  some  arc  is  unessential,  it  remains  unessential  through  the 
end  of  the  computation.  An  unessential  arc  cannot  belong  to  a  maximum  surplus  cut.  We 
show- that  few  iterations  are  enough  to  make  a  new  arc  imessential. 

Lemma  5  From  the  iteration  i  +  2  on,  the  (current)  capacity  of  cut  {Si,Ti)  is  greater  than 
Si+i. 

Proof.  At  the  beginning  of  iteration  t  + 1  the  capacity  of  cut  (Sj,  T*)  is  equal,  by  definition, 
to  Hi,  the  remaining  total  demand.  In  iteration  t  + 1  the  capacity  of  each  cut  decreases  by 
at  most  Hi  —  Hi+i  and  then  increases  by  (5,-+i/Si+i)P,  where  B  is  the  weight  of  this  cut. 
It  means  that  in  iteration  t  +  1  the  capacity  of  cut  {Si,Ti)  first  decreases  but  does  not  go 
below  Hi^\.  Then  it  increases  by  more  than  Hi+i,  because  its  weight  is  greater  than  Bi+i- 
Therefore  at  the  end  of  iteration  t  +  1  the  capacity  of  {Si,Ti)  is  greater  than  2Hi+i.  From 
now  on  the  capacity  of  (5*,ri)  is  always  greater  than  5i+i,  because  the  total  decrease  of 
this  capacity  cannot  be  greater  than  the  total  demand  left  after  iteration  t  +  1,  whidi  is 

Hi+i- 

Putting  the  above  argument  in  a  more  formal  way  we  get  the  following  bound  on  the 
capacity  of  {Si,  Ti)  at  the  beginning  of  iteration  t+l,  where  I  >  2.  According  to  our  notation, 
the  capacity  function  at  the  beginning  of  iteration  t  +  /  is  • 


> 


rO  +  {6i+i  -  6i)Bi  -{Hi-  Hi+i.i )  . 

{6i+i  —  6i+i)Bi  +  (5i+i  —  6i)Bi  —  Hi  +  Hi+i-i 

{6i+i  —  6i+i)Bi  +  Hi+i-i 

(^i+2 
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>  (^*+2  — 

=  -ETi+i  • 


We  first  prove  an  O(mlogm)  bound  to  show  the  mzdn  idea.  To  prove  an  0{m)  bound, 
we  will  need  a  finer  accounting  strategy. 

Theorem  6  Newton’s  method  solves  MCVC  in  0{mlogm)  iterations. 

Proof.  We  use  Inequality  (4)  to  show  that  after  O(logm)  iterations  a  new  arc  becomes 
imessentisd.  Let  the  current  iteration  be  the  zth  one.  Let  I  =  [log  mj  +  2.  It  follows  from 
Inequality  (4)  that 

Ei+iBi+i  <  — (10) 

771* 

If  Si+i  <  ^Bi,  then  there  exists  an  arc  e  e  (Si.Ti)  such  that  6(e)  >  Bi+i-  Such  an  arc 
is  unessential  from  iteration  i  +  I  on. 

If  Bi+i  >  ^Bi,  then  Bi+i  >  as  well,  and  Inequality  (10)  implies  that  Ei+i  < 

^Ei+i.  Lemma  5  says  that  the  capacity  of  (5*,  T*)  at  the  beginning  of  iteration  i  +  /  +  1 
is  greater  than  Ei+i,  so  its  greater  than  mEi+i-  It  means  that  the  capacity  of  some  arc  in 
{Si,Ti)  in  (which  is  the  network  at  the  beginning  of  iteration  t  +  Z  +  1)  is  greater 

Ei+i,  the  total  demand  in  ■  Such  an  arc  is  unessential  from  iteration  t  +  Z  +  1 

on.  ■ 


Theorem  7  Newton’s  method  solves  MCVC  in  0{m)  iterations. 


Proof.  Lemma  1  implies  that  for  each  iteration  t, 

Bi 


or 


Bi.i  -  2 


Ei.i  -  2 


(11) 


(12) 


We  first  bound  the  number  of  iterations  i  for  which  (11)  holds.  Let  (ftt)f-i  be  the  sequence 
of  jBi’s  in  these  iterations  and  be  the  positive  fc^’s  arranged  in  nonincreasing  order. 

It  is  easy  to  see  that  sequences  and  satisfy  the  conditions  of  Lemma  6,  so 

q  <p<rn. 


Now  we  bound  the  number  of  iterations  i  for  which  (12)  holds.  Here  the  argument  is 
more  involved,  because  Hi’s  are  not  just  subsums  of  a  set  of  tti  elements.  They  are  related 
to  the  current  demands  and  capacities,  which  vary  from  iteration  to  iteration- 
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To  avoid  towering  subscripts,  we  reniunber  iterations  taking  into  an  account  only  itera¬ 
tions  with  (12).  It  means  that  now  iteration  i  is  what  used  to  be  the  ith  iteration  with  (12). 
AB  subscripts  refer  to  the  new  numbering. 

We  assign  iterations  to  arcs  in  such  a  way  that  at  most  5  iterations  are  assigned  to  one 
arc.  We  stop  the  process  of  assigning  when  all  but  at  most  g  -f  2  =  [log  m]  +  5  iterations 
have  been  assigned.  Assume  that  there  are  still  at  least  5  ■+  3  unassigned  iterations.  Let 
iteration  i  be  the  first  among  them.  Consider  cut  [Si,Ti).  Let  p  be  its  cardinality.  Let  for 
1  <  j  <  p  and  l<l<q,  'fij  be  the  capacity  of  the  jth  arc  in  {Si,Ti)  (assuming  any  order 
on  the  arcs  in  (5*,  Ti))  at  the  beginning  of  iteration  i  +  l  +  2.1t  follows  from  Lemma  5  that 
for  each  !<!<?> 

(13) 

For  each  arc,  the  difference  between  its  capacities  at  the  beginning  of  iterations  1:  +  1  and 
lb  +  2  is  not  greater  than  the  total  demand  at  the  beginning  of  iteration  fc  -f  1,  which  is  Ek- 
It  means  that 

< -S’x-i-J-i-i  <  (14) 

The  definition  of  q,  (13),  and  (14)  imply  that  matrix  (71  j/iTi+i)  satisfies  the  conditions  of 
Lemma  7.  It  means  that  there  exists  lb  >  3  such  that  at  least  k -2  elements  from 7^,1,  7it,2, 
7fc,p  are  greater  than  l/2*'5^i+i  >  Ei+k+i-  Therefore  the  arcs  corresponding  to  these 
lb  —  2  elements  are  unessential  at  iteration  i  -1-  fc  -h  2.  We  assign  iterations  i  through  i  +  k  +  1 
to  these  arcs.  Notice  that  none  of  these  arcs  was  unessential  at  iteration  i,  so  none  of  the 
previous  iteration  was  assign  to  any  of  them. 

During  this  process  of  assigning  iterations  to  arcs,  no  more  than  5  iterations  are  assign 
to  one  arc,  and  no  more  than  O(logm)  iterations  are  left  unassigned.  We  conclude  that 
there  are  at  most  5m  +  O(logm)  iterations  for  which  (12)  holds.  ■ 

Lemma  6  Let  cci  >  Q2  >  •  •  •  >  Op  >  0  and  fii  >  fi2  >  -  ••>  ftq  >  0  be  such  that 
1-  /it-fi  <  |/*t.  for  i  =  1,2, . .  .,g  -  1, 

2.  Hg  >  Op, 

3.  Pi  <  ^  Pt}.fort  =  1,2,. ..,g. 

Then  g  <  p. 
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Proof.  Let  &j  =  aj  +  a^+i  H - 1-  Op.  Condition  3  implies  that  n\  <  a\.  Thus,  to  prove 

g  <  p,  it  is  sufficient  to  show  that  each  of  the  intervals  (0,5:^],  (op.Qp-i],  (03,02], 
(q2,  oi]  contains  at  most  one  element  from  (/i,-).  Condition  2  implies  that  only  the  last 
element  can  be  in  (0,  Opj.  Let  for  some  1  <  J  <  P  —  1  aJid  —  1,  A*t€  (oj+i,  Oj). 

We  show  that  pt+i  <  *j+i- 

If  Q,+i  >  jOj,  then 
If  Oj+i  <  |qj,  then 

W+i  ^  -  \^j  -  ~  ^ 

The  above  inequality  and  Condition  3  imply  ^+1  <  ^j+i  ■  ■ 

Lernina  7  Let  (oij)  be  a  g  x  p  matrix  such  that 


1-  ?  >  log?  +  3, 

2.  the  sum  of  each  row  is  not  less  than  1, 

3.  |oi+ij  -  Oijl  <  l/2\  for  1  <  t  <  g  and  1  <  j  <  p. 


If  Oij  >  1/2*,  then  we  call  this  and  all  subsequent  elements  in  column  j  good  elements.  There 
is  ife  >  3  such  that  row  k  contains  at  least  A  —  2  good  elements. 


Proof.  Condition  3  implies  that  if  1  <  t'  <  t"  <  g  and  1  <  j  <  p,  then 

1 

ff  Oij  is  the  &st  good  element  in  column  j  and  i  >  2,  then 

1  J_  ,  _L.  -  J_ 

Otj  £  J  +  2»-i  ~  2*~^  2*“'  2*”^  * 


(15) 


(16) 


Assume  that  for  all  A  >  3,  row  A  contains  at  most  A  —  3  good  elements.  We  will  get 
contradiction  by  showing  that  the  sum  of  the  last  row  is  less  than  1.  Let  ii,i2>  •  •  •  j  J/ 
the  indices  of  all  columns  with  at  least  one  good  element.  Let 

the  first  good  elements  in  these  columns.  Let  •  •  •  »il  l>e  ordered  in  such  a  way  that 

*1  <  *2  <  - '  •  <  Bow  ik  contains  at  least  A  good  elements  (ott  Ji » » — » ^kJt ) 

ifc  >  A  +  3.  (^'^) 
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The  Sinn  of  the  last  row  is  at  most 


29 


+  “9.5i  +  <^0.15  + - !■  O:. 


5.J2 


*■<1,31 


<  ^  +  (“ij  +  jiT^)  + 

1,1  1  , 

^  8  +  (^  +  +  • 


<  1. 


■  +  (a, 


*iji 


,  1 

+  (2iP2 


+ 

+ 


-r— ) 

2*f-i  ’ 


The  first  inequality  follows  from  (15),  the  second  one  from  (16)  and  Condition  1,  and  the 
third  one  from  (17).  ■ 


5  Concluding  Remarks 

We  proved  that  Newton’s  method  solves  a  linear  fractional  combinatorial  optimization  prob¬ 
lem  in  a  strongly  polynomial  number  of  iterations.  We  did  it  by  combining  the  fast  con¬ 
vergence  of  the  method  with  the  limiting  way  the  numbers  involved  in  an  LFCO  problem 
can  be  spread  along  the  real  line.  That  “limiting  spread”  follows  from  the  fact  that  these 
numbers  are  obtained  using  only  few  multiphcations/divisions.  Observe  that  for  example 
the  generalized  flow  problem  does  not  have  this  property,  and  we  do  not  know  how  our  t3rpe 
of  analysis  could  be  useful  in  such  a  case. 

We  showed  that  Newton’s  method  gives  an  0(m)  iteration  algorithm  for  the  transship¬ 
ment  problem  when  the  maximum  arc  cost  is  minimized.  One  iteration  in  this  algorithm  is 
dominated  by  the  maximum  flow  computation,  so  the  overall  running  time  is  0{rr?n).  This 
improves  the  upper  bound  achieved  with  Megiddo’s  parametric  search.  We  believe  that  the 
bound  on  the  number  of  iterations  is  0(n),  but  we  have  not  been  able  to  prove  it. 

We  may  have  a  transshipment  network  with  two  cost  functions,  one  related  to  minimizing 
the  sum  of  arc  costs  and  the  other  to  TniwiTniwTig  the  maximum  arc  cost .  In  the  same  0  (m^n) 
time  bound  we  can  perform  the  “double  minimization”.  To  minimize  the  sum  of  the  arc 
costs  over  the  transshipmmts  which  minimize  the  maximum  arc  cost,  we  first  minimize  the 
maximum  arc  cost  and  get  in  this  way  additional  capacity  constraints.  Then  we  minimize 
the  sum  of  arc  costs  taking  into  an  account  these  additional  constraints.  To  minimize  the 
maximum  arc  cost  over  the  transshipments  whidi  minimize  the  sum  of  arc  costs,  we  first 
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compute  a  transshipment  which  minimizes  the  sum  of  arc  costs,  and  its  dual,  and  identify 
all  arcs  which  have  to  be  saturated  in  all  such  transshipments.  Then  we  can  minimize  the 
maTiTmim  arc  cost  of  the  remaining  axes. 
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